Empirical tight binding parameters for GaAs and 
MgO with explicit basis through DFT mapping 
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Abstract. The Empirical Tight Binding(ETB) method is widely used in atomistic 
device simulations. The reliability of such simulations depends very strongly on the 
choice of basis sets and the ETB parameters. The traditional way of obtaining the 
ETB parameters is by fitting to experiment data, or critical theoretical bandedges and 
symmetries rather than a foundational mapping. A further shortcoming of traditional 
ETB is the lack of an explicit basis. In this work, a DFT mapping process which 
Q . constructs TB parameters and explicit basis from DFT calculations is developed. The 

method is applied to two materials: GaAs and MgO. Compared with the existing 
TB parameters, the GaAs parameters by DFT mapping show better agreement with 
the DFT results in bulk band structure calculations and lead to different indirect 
valleys when applied to nanowire calculations. The MgO TB parameters and TB basis 
■ functions are also obtained through the DFT mapping process. 
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Modern semiconductor nanodevices have reached critical device dimensions in the 
range of several nanometers. These devices consist of complicated two and three 
dimensional geometries composed of multiple materials. Typically, about 10000 to 
10 million atoms are in the active device region with contacts controlling the current 
injection. This finite extent suggests an atomistic, local and orbital-based electronic 
structure representation. Quantitative device design requires the reliable prediction of 
bandgaps and band offsets within a few meV and effective masses at principal symmetry 
points within a few percent. 

Ab-initio methods that have no adjustable parameters offer such atomistic 
representations. However, accurate models such as hybrid functionals [1], GW [2] 
and BSE approximations [3] are computationally far too expensive to be applied on 
multi-million atom devices. More approximate ab-initio methods such as the local 
density approximations (LDA) and generalized gradient approximations (GGA) [4] do 
not reproduce band gaps, relative band offsets, and effective masses accurately enough. 
Empirical methods such as the empirical tight binding (ETB) method are numerically 
much more efficient. The accuracy of ETB is hereby limited by the parameters fitting. 
Previous ETB simulations in semiconductor nanodevices such as resonant tunneling 
diodes [5], quantum dots [6] and strained Si/SiGe quantum wells [7] showed quantitative 
agreement with experiments. 

The accuracy of the ETB method depends critically on the careful calibration 
of the empirical parameters. The common way to determine the ETB parameters 
is to fit ETB results to experimental band structures [8] [9]. One shortcoming of 
this method is its requirement of experimental data that are often not available for 
new and exotic materials. In addition, the ETB basis functions remain unknown, 
which makes it notoriously difficult to predict wave function dependent quantities with 
high precission. To overcome these shortcomings, some approaches were developed 
to construct localized basis functions from ab-initio results such as localized wannier 
functions [10] or quasi-atomic orbitals [11] [12] [13]. Unfortunately, these functions are 
either not reliably centered at atoms, or resulting Hamiltonian requires long distance 
coupling with large number of neighbors, which is numerically expensive. Nanoelectronic 
devices are increasingly based on complex heterstructures in 1 or 2 dimensions including 
atomistc disorder and strain. In this realm models with one or two neighbors are 
conceptually preferrable and simple to implement. In contrast the validity of many 
neighbor bulk basis states are very questionable in this domain. 

In this work, a DFT mapping method that constructs ETB parameters from ab- 
initio calculations is presented. This method allows to determine ETB basis functions 
that are centered at atoms and it limits the interatomic coupling to the first or second 
nearest neighbors. Since the method does not require experimental results, the ETB 
parameters are less empirical. Two materials are considered in this work: 1) the well 
known GaAs to validate the method and 2) MgO which is recently used as a magnetic 
tunneling barrier material in Magnetoresistive Random- Access-Memory devices. MgO 
lacks the elaborate experimental analysis, but it is known to have small spin-orbit 



3 



interaction and a large band gap [14]. 
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Step 3: Low rank approximate H DFT to H TB 
Approximate H TB by two center integrals 
Omitting coupling > 1st NN 
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Step 6: Get two center integrals from step 3 
and basis functions from step 2 





Figure 1. The process of TB parameters construction from DFT calculations. 



Figure 1 shows the flow chart of the mapping method. The first step is to perform 
ab-initio calculations of the band structure of a material. In general, any method that 
is capable to calculate electronic band structures and wave functions is suitable here. 
However, in this work, DFT calculations including hybrid functionals corrections for 60 
electronic bands performed with VASP version 5.2 [15] are used. Inthe second step, the 
ETB basis functions for each type of atom are defined as 

# n>J>Tn (r) = (r, 9, 4>) = (r) Y hm {9, 0) , (1) 

where the functions Y\ %m are the complex spherical harmonics and the functions R Ui i are 
exponentially damped plane waves 

N 

Rn,l ( r ) = ^ Sil1 ( Air ) +bi COS (^ r )] T ' n ~ 1 eXP ( _ttir ) • ( 2 ) 
i=l 

Here, r = (r, 9, (p) is the position vector with the respective atom centered at the origin, 
I, m are the angular and magnetic quantum numbers of the orbital basis function, 
and n is the principal number of the atomic orbital as used in Hiickel type basis 
functions [16]. The remaining parameters Oj, &i, (Xi, Aj are the fitting parameters. The 
ETB basis functions are spin independent. With a given set of ETB basis functions 

in the third step, a transformation matrix between these chosen DFT basis functions 
^DFT anc ^ ^ e ^TB ^ s c& l ciua t e d- Since the number of the ETB basis functions is 
smaller than the DFT basis funcitons, this transformation matrix is rectangular and 
represents a low rank approximation. [17] Then, the DFT Hamiltonian is transformed to 
the tight binding representation. The ETB Hamilton matrix elements are approximated 
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by two center integrals according to the Slater-Koster table [18] [19]. Any non-zero off- 
diagonal element of the overlap matrix is neglected. ETB Hamilton matrix elements 
beyond either 1st or 2nd nearest neighbor coupling are neglected. In Step 4, the band 
edges, effective masses and eigen functions of the Hamiltonian at high symmetry points 
are calculated and compared to the corresponding DFT results. The overlaps of the 
ETB basis functions are also determined. In the fifth step, all fitting parameters are 
adjusted to improve 1) the agreement of the ETB results with the DFT ones and 2) 
reduce the overlap matrix of the ETB basis functions to the unity matrix. Steps 2-5 
are repeated until the convergence criterion is met, i.e. when the maximum difference of 
DFT and ETB band edges are within 10 meV, when the effective masses agree within 
5% and when the eigenfunctions of DFT and ETB calculations agree by at least 90%. 
Step 6 requires to extract the converged ETB basis functions and the ETB two center 
integrals. 
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Figure 2. Band structure and density of states of GaAs by DFT(green dashed lines), 
TB using parameters in ref[9] and TB using parameters by this work. 

For GaAs and MgO, eigenfunctions and eigen energies of the lowest 16 bands of the 
DFT calculations from L to T and r to X are taken into account for the ETB fitting 
method. Here, the wave functions for the topmost valence bands and lowest conduction 
band valley are considered in the fitting of ETB eigenfunctions to the DFT ones. The 
overlap of the ETB basis functions is partly minimized in the fitting process. Most of the 
overlap matrix elements vanish, but the maximum overlap i.e. in this case the overlap 
of the p orbitals of cations and d orbitals of anions, e.g. {Px a \dy s z ) = 0.86 at the T point 
remains comparably high. GaAs is parameterized for the 1st nearest neighbor sp 3 d 5 s* 
ETB model. The resulting parameters can be found in table 2. The band structure and 
density of states of GaAs are shown in Fig. 2. Calculated bandstructures of the DFT 
method, of the ETB method with parameters of the present mapping method and of 
the ETB method with previously published parameters [9] are compared in Fig. 2 as 
well. The ETB calculations with new parameters agree very well with the DFT results 
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Figure 3. Comparison of E-k diagram of conduction bands (left figure) and 
Transmission (right figure) of GaAs nanowire . The red lines are results using TB 
parameters generated by this work; the blue line are results using TB parameters by 
ref[9]. 









GaAs 






MgO 




Quantity 


DFT 


TB 


Error 


ref[9] 


DFT 


TB 


Error 


E 9 (r) 


1.420 


1.449 


2.96% 


1.413 


7.831 


7.499 


4.2% 


E g {X) 


1.973 


1.947 


0.9% 


1.898 


12.161 


11.819 


2.8% 


E g (L) 


1.728 


1.718 


0.6% 


1.714 


10.871 


10.469 


3.7% 




0.0692 


0.0737 


6.5% 


0.0657 


0.396 


0.458 


15.6% 


m *x,i) 


1.140 


1.117 


2.0% 


1.881 








m *x,t) 


0.219 


0.231 


5.5% 


0.175 








m (L,l) 


1.700 


1.756 


3.3% 


1.728 








m (L,t) 


0.133 


0.138 


3.8% 


0.098 








7i 


6.964 


6.985 


1.1% 


7.388 


0.952 


0.889 


6.6% 


72 


2.084 


2.151 


3.6% 


2.367 


0.277 


0.219 


20.9% 


73 


2.972 


2.980 


1.1% 


3.098 


0.376 


0.234 


37.7% 



Table 1. Comparison of important bandedges and effective masses of GaAs and MgO 
by DFT and TB. 

for energies below 5 eV whereas ETB calculations with previously published parameters 
deviate from the DFT results already at about 2 eV. Relevant band edges and effective 
masses are compared in table 1, demonstrating a much better fit. 

Figure 3 compares results of ETB calculations of the band structure and the 
transmission coefficient of a GaAs 3.5 nm x 3.5 nm squared nanowire when the new and 
the previously published parameter set are used. The results agree for the T point at 
lower energies, but they deviate significantly for the indirect conduction valley at about 
1.9 eV. The difference in the confinement energy of this conduction valley originates 
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Figure 5. Radial part of Basis functions of Mg and atoms in MgO. 



from different transverse effective masses at the L point of the two ETB parameter sets 
(see Table 1). The modeling go the details of the high symmetry points is particularly 
important for the emerging concept of gamma-L transistors [20]. 

MgO cyrstalizes in Rocksalt structure. Each Oxygen atom has 6 Magnesium atoms 
as 1st nearest neighbors and 12 Oxygen atoms as 2nd nearest neighbors. The valence 
bands of MgO are formed by hybridized orbitals of Oxygen atoms: The s-orbitals and 
]?-orbitals of Oxygen contribute to valence bands around —17 eV and — 1 eV respectively, 
while the s*, (i-orbitals of Oxygen and orbitals of Magnesium contribute to the band 
structure for energies beyond 7 eV (see Fig. 4). According to this, MgO is parameterized 
for a 2nd nearest neighbor sp 3 d 5 s* ETB model (for parameters see table 3). Within this 
model, the interaction between two Oxygen atoms is required to produce the correct 
valence bands. The interaction between two Magnesium atoms is omitted, since the 
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onsite energies of Magnesium orbitals are higher than 20 eV so that the omission will 
only affect bands close to that high and technically irrelevant energy. The ETB band 
structure matches the DFT result well within the energies 5 to 20 eV. Important band 
properties calculated in the DFT and the ETB method are listed in Table 1. The 
basis functions of Oxygen and Magnesium are shown in Fig. 5. The basis functions 
of Oxygen are more localized while the basis functions of Magnesium are more plane 
wave like functions. This difference in the localization originates from the fact that 
Magnesium orbital energies are higher than the Oxygen ones. 

In conclusion, a method to determine ETB parameters from density function 
theory calculations is developed. The method is applied to GaAs and MgO. First 
nearest neighbor ETB parameters and basis functions for an sp 3 d 5 s* model of GaAs are 
presented. Results of this parameterization agree well with the DFT calculations. The 
new ETB parameters lead to lower indirect conduction bands when applied to GaAs 
nanowires. Second nearest neighbor ETB parameters and basis functions for an sp 3 d 5 s* 
model of MgO are also obtained. The ETB results with this parameterization also agree 
well with the DFT calculations. 

nanoHUB.org computational resources operated by the Network for Computational 
Nanotechnology funded by NSF are utilized in this work. The research was funded by 
the Lockheed Martin Corporation and NSF (Award No. 1125017) 

Table 2. Tight Binding Parameters for bulk GaAs 



Parameter 


value 


Parameter 


value 


E sa 


-4.5863 


S*aPc° 


2.6877 


E 


1.4694 


S* c PaO 


1.8335 


E 

J ->s*a 


10.0480 


s a d c a 


-2.1172 


Eda 


11.2878 


s c d a a 


-2.9128 


E sc 


-1.3323 


s* a d c a 


-0.4974 


E 

J -<pc 


9.5885 


s* c d a a 


-2.9971 


E 


25.6752 


ppa 


3.8065 


Edc 


35.2863 


ppn 


-1.5010 


A a 


0.1259 


Pad c o 


-1.2077 


A c 


0.1235 


Pcd a cr 


-1.9855 


ssa 


-1.7615 


p a d c 7T 


3.1547 


s*s*a 


-0.8374 


Pcd a ii 


2.3234 


s* a s c a 


-1.1173 


dda 


-1.9986 


s a s*a 


-2.9313 


ddn 


3.1681 


SaPcO- 


2.1768 


dd5 


-2.3137 


ScPaCT 


3.6705 
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Table 3. Tight Binding Parameters fo r bulk MgO 



Parameter 


Value 


Parameter 


Value 


E sa 


-7.1496 


Pad c ir 


1.5284 


E 


5.8926 


PcdaK 


-4.0453 


E 

J -'s*a 


23.8138 


d a d c a 


-1.0038 


Eda 


40.0285 


d a d c n 


5.0830 


E sc 


38.4754 


d a d c 5 


-0.6323 


E 


32.1465 




-0.2718 


F, 


45.5084 




-0.4690 


Edc 


57.9865 


s a s* a (J 


-0.0001 


K 


0.0031 


SaPaO 


0.3388 


A c 


0.0298 


S*aPa(? 


0.1965 


s a s c a 


-0.1192 


s a d a a 


-0.3380 




1.6477 


s* a d a a 


-0.4407 




-0.6008 


PaPaO- 


0.4371 


s a s*a 


-0.6347 


PaPaK 


-0.0641 


SaPcCT 


-2.0013 


PadaCT 


-0.4039 


ScPaO- 


0.3283 


PadalT 


0.6986 


S* a Pc<7 


2.1584 


d a d a a 


-2.3768 


S* c Va° 


2.1282 


d a d a n 


0.4556 


s a d c a 


0.6641 


d a d a 5 


0.0967 


s c d a a 


-2.9483 






s* a d c a 


1.6890 






s* c d a a 


3.1534 






PaPc^ 


0.1743 






PaP C K 


-0.4703 






Pad c cr 


-1.9960 






p c d a a 


0.0519 
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